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ABSTRACT 


Shock structure during ionization of a hydrogen- helium mixture has been 
followed using hydrogen line and continuum emission measurements. A reac- 
tion scheme is proposed which includes hydrogen dissociation and a two-step 
excitation-ionization mechanism for hydrogen ionization by atom-atom and 
atom-electron collisions. Agreement has been achieved between numerical 
calculations and measurements of emission intensity as a function of time for 
shock velocities from 13 to 20 km/sec in a 0, 208 H 2 - 0. 792 He mixture. The 
electron temperature was found to be significantly different from the heavy 
particle temperature during much of the ionization process. Similar time- 
histories for H(3 and continuum emission indicate upper level populations of 
hydrogen in equilibrium with the electron concentration during the relaxation 
process. The expression for the rate constant for excitation of hydrogen by 
atom-atom collisions that best fit the data was 

k AA = 4 -° x 10 " 17 (^f) l/2 ex P (-1%*) cm 3 sec* 1 , 

where it has been assumed that the excitation cross section is the same for 
hydrogen and helium collision partners. The electron-atom excitation rate 
constant. 


k 3 = 7. 5 X 10 


“CSf 


exp (- 10 /kT > cm^ sec 


determined from this investigation, was in agreement with recent electron 
beam cross section measurements. 
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I. INTRODUCTION 



While shock structure during ionization 
relaxation has been a subject of recent research 
for air, the noble gases, and akali-metal vapors, 
the knowledge of hydrogen ionization relaxation is 
limited. An understanding of shock structure in 
hydrogen- helium mixtures is needed to determine 
the importance of nonequilibrium flow during 
atmospheric entry of probes to the outer planets. 
The presence of a thick nonequilibrium layer in 
the entry body flow-field could result in a signifi- 
cant reduction in the estimated radiative heating 
tnat the probe must be designed to withstand. In 
the present study, shock tube measurements and 
numerical analysts have been performed in order 
to understand the ionization and thermal processes 
which govern the shock structure in hydrogen- 
helium mixtures. 

Much of the knowledge of collision processes 
governing the ionization of monatomic gases has 
come from studies of argon ionization by Petschek 
and Byron (Ref. 1), Harwell and Jahn (Rei. 2), 
Kelly (Ref. 3), Wong and Bershader (Ref. 4), 
Hoffert and Lien (Ref. 5), and others. Argon 
ionization is initiated by the excitation of argon 
atoms by atom-atom collisions. The excited 
atoms are then quickly ionized due to subsequent 
collisions. After a minimum number of electrons 
have been produced Sy atom-atom processes, the 
ionization is governed by a two-step electron-atom 
collision process. Good agreement has been ob- 
tained between experimental data and calculations 
based on this reaction scheme, assuming the 
existence of a separate electron and heavy par- 
ticle temperature. The electron temperature was 


found to be coupled to the ionization relaxation 
because ionization by electron-atom collisions was 
the dominant form of electron energy transfer. 

The ionization and thermal processer govern- 
ing shock structure in hydrogen-helium mixtures 
have been postulated (Ref. 6) based on the results 
of the investigations of argon. A reaction scheme 
is proposed which includes hydrogen dissociation, 
electronic excitation, and ionization by atom- atom 
and atom-electron collisions. Numerical solu- 
tions are obtained for the variation of fluid vari- 
ables and species concentrations behind strong 
shock waves where the electron temperature has 
been allowed to differ from the heavy particle 
temperature during the relaxation process. 

Ionization relaxation has been followed using 
spectroscopic emission measurements in a conical 
driver electric arc shock tube. Hydrogen line and 
continuum emission intensities were obtained as a 
function of time behind the incident shock wave in 
a 0. 2079 H£ - 0. 7921 He mixture. This mixture 
allowed temperatures of interest .o be achieved 
within the performance limits of the shock tube. 
From the agreement of calculated emission inten- 
sity time histories with the experimental data, it 
was inferred that a model describing the important 
processes governing shock structure had been ob- 
tained. New estimates of the rate constant for the 
excitation of hydrogen by atom-atom collisions 
were obtained and the rate constant for excitation 
of hydrogen by electron-atom collisions was found 
to be in agreement with electron beam cross sec- 
tion measurements. 
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II. EXPERIMENT 


The time history of hydrogen line and con- 
tinuum emission has been measured behind the 
incident shock wave in a mixture of hydrogen and 
helium. The effect on radiation relaxation rate of 
variations in shock velocity between 13 and 
20 km/sec and of initial pressures between 33. 3 
and 267 N/m^ (C. 25 and 2.0 torr) has been deter- 
mined. Shock waves were studied in a 0. 152-m- 
diameter shock tube using the recently developed 
conical electric arc driver described by Menard 
(Ref. 7). Helium driver gas behind a thin alumi- 
num diaphragm Is heated and accelerated by the 
arc discharge of a 290 - kilojoule capacitor bank. 

Emission measurements were made of the 
Hp line and the hydrogen continuum using a Jar -ell 
Ash, 0. 75-m, grating spectrometer. The spec- 
tral interval of the continuum channel was selected 
to avoid contributions from strong driver gas 
radiators. Channels were also selected to observe 
radiation from a helium line, a portion of a hydro- 
gen molecular band, and a line of chromium which 
is a major radiator in the driver gas. The spec- 
tral intervals of each channel and a schematic 
drawing of the instrumentation are shown in Fig. 1. 
The photomultipliers (RCA 1P28A) and electronics 
had a combined response time of less than 0. 2 
fisec. The time resolution of the optical system 
was less than the photomultiplier response time 
for the shock velocities studied in this investiga- 
tion. 


Arrival of the shock wave at the spectrometer 
station was detected with a Kistler pressure 
transducer, and shock velocities were measured 
using a series of photomultipliers and thin film 
resistance gages. The test gas used was a 
Matheson Gas mixture of 20. 79% hydrogen with 
the balance helium. The mixture was prepared 
using ultra-high-purity grade gases and analyzed 
to have less than 2 ppm water, 0. 5 ppm hydro- 
carbons, and 6 ppm nitrogen. Before each run the 
shock tube was evacuated to a pressure of at least 
1. 33 X 1C 3 N/m 2 (1 X 10*5 mm H g) The appar- 
ent leak rate was typically 4. 4 X 10" 4 N/m 2 -sec 
(2 X 10“ 4 torr/min). After evacuation the test 


gas was continuously flushed through the shock 
tube for 2 min at the desired fill pressure, 

Shock tube runs were made over the range of 
shock velocities from 13 to 20 km/sec and for 
initial pressures of 33, 133, and 267 N/m^ (0, 2$, 
1.0, and 2,0 torr). A sample of the data obtained 
from a run with 15.0 km/sec shock velocity Is 
shown in Fig. 2. The slope of the signals from 
both the Hp and hydrogen continuum channel (\ = 
514. 5 mn) were zero immediately behind the shock 
wave, increased to a maximum value, and then 
returned to zero at. the steady state signal levels 
were approached. A characteristic relaxation 
time for the continuum signal, t^, was defined 
as shown in Fig. 2 by the intersection of the max- 
imum slope line and the steady state value. The 
relaxation time for the Hp emission was equal to 
the continuum relaxation time within the measur- 
ing accuracy. Both tj-jj and the time between 
shock arrival and the initial jump in the intensity 
signal (induction time) decreased with increasing 
shock velocity, although the induction time became 
a smaller fraction of t^ &t the larger shock 
velocities. For shock velocities greater than 
17 km/sec, the Kistler g*< y e was not able to re- 
solve the arrival of the shock wave and the initial 
jump in the intensity signal was assumed to indi- 
cate shock arrivals 

From the relative intensities and time varia- 
tions of the He and band channels it was con- 
cluded that the hydrogen continuum was the domi- 
nant source of radiation In these spectral inter- 
vals. Thus the radiation in ail channels observed 
was governed by the hydrogen ionization process 
with no detectable sign of H ^ dissociation. The 
signals from the chromium channel were similar 
to the continuum channel until the arrival of the 
contact surface, when the signal level would in- 
crease by a large factor. The sudden decrease in 
the Hp intensity at the same time as the Jump in 
the Cr channel indicates that a test slug free of 
driver gas contamination was formed in the shock 
tube. Test times \yere more than adequate for 
-relaxation measurements over the range of condU 
tions used, and slug lengths were found to be 
approximately the tube diameter. 
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in. DESCRIPTION OF ANALYSIS 


i 



A, Reaction Scheme and Rate Equations 

The reaction scheme describing the important 
colHsionfii processes in hydrogen-helium ionizing 
shock waves has been modeled after the results 
obtained from studies of argon ionization (Refs. 
4-**), Reactions are included which dissociate the 
molecular hydrogen, excite the electronic states 
of hydrogen and helium, and ionize the hydrogen 
and helium by collisions with atoms and with elec- 
trons. Table 1 lists the 11 steps of the reaction 
scheme. Photo -ionization and precursor effects 
are considered negligible, consistent with the 
findings of Bates, Kingston, and McWhirter 
(Refs. 8 and 9). 

Ionization of hydrogen is initiated by atom- 
atom Collisions which produce atoms in electronic- 
ally excited states that are subsequently ionized 
by additional collisions. A simplified model is 
assumed in which hydrogen in the first excited 
state is produced by collisions with hydrogen 
(reaction 5) and from collisions with helium (reac- 
tion 6). It is assumed that excitation to the n = 2 
level is the rate -limiting step and the excited 
atoms are rapidly ionized by further collisions. 
Hydrogen ionization by electron-atom collisions 
also proceeds by a two-step excitation- ionization 
process in which the ionization of excited hydrogen 
is assumed to be rapid. The direct ionization of 
hydrogen by electron-atom collisions is included 
as reaction 1. Helium is ionized by collisions 
with electrons from a single step process in reac- 
tion 2 and by an excitation-ionization process in 
reaction 4. 

The reaction scheme for the dissociation of 
hydrogen has been obtained from shock tube inves- 
tigations (Refs. 10-14). Hydrogen is first dis- 
sociated bv collisions with itself and helium; but, 
subsequently, it is dissociated by collisions with 
atomic hydrogen, ions, and electrons. 

The knowledge of rate constants for hydrogen 
ionization is incomplete. No definitive estimates 
of atom-atom rate constants have been made. 

The rate constants for excitation and ionization by 
electron-atom collisions were calculated using 
data obtained from electron beam cross section 
measurements. A fit of the cross section data 
was obtained by matching the initial slopes of the 
cross section data and the expression = 
ai (1 - Ei q/E), Assuming a Maxwell- Boltzmann 
distribution of reactants, the resulting expression 
for the rate constants (Ref. 1 5) Is 

where E t q is the threshold energy for reaction 
(&), is ’an adjustable parameter, and p is the 
reduced mass. Initial estimates of the rate con- 
stant for atom-atom excitation of hydrogen were 
made assuming that the ratio of the atom-atom 
and atom-electron excitation cross section for 
hydrogen was approximately the same as the ratio 
for the Noble gases (Ref. 16). The rate constant 
expressions used in this investigation are shown 
in Table 2. 
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Hydrogen dissociation reactions have been 
studied in shock-heated mixtures of argon and 
hydrogen over the range of temperatures 
5000 > T > 2500° by Myerson and Watt (Ref. 10), 
Jacobs, Giedt, and Cohen (Ref. 11), Patch (Ref. 
12), Rink (Ref. 13), and Sutton (Ref. 14). The 
expressions for the rate constants given by Jacobs 
et al. were found to represent the best compro- 
mise of the experimental data obtained by all 
investigators and were suitable for extrapolation 
to temperatures greater than 5000 °K. The rate 
constants for the hydrogen dissociation reaction 
with helium and with argon as collision partners 
were assumed to be equal. 

The rate constants for the reverse direction 
of each reaction step were found, assuming the 
principle of detailed balancing. Thus, the reverse 
rate constant for a given step was equal to the 
rate constant for the forward direction divided by 
the equilibrium constant for that reaction. The 
equilibrium constants were obtained by fitting 
data produced by the thermochemistry program 
of Menard and Horton (Ref. 22) with the expres- 
sions 

K h + = 4. 05 X 10* 9 T 3 / 2 

X exp (-1. 578 X 10 Vt) mole/cm 3 , 

K„ + = 1. 62 X 10' 8 T 3 / 2 

He 

X exp (-2. 853 X n S /T) mole /cm 3 , 

and 

K h = 3. 7 [l . 0 - exp (-1. 50 X 10 8 T" 2 )] 

X exp mole/cm 3 . 

Six equations are needed to describe the con- 
centrations of the species H 2 , H, H*, He, He + , 
and e in the present system. Reaction rate equa- 
tions for H, e, and He + are used along with 
three conservation of atomic species relations. 

The concentrations of all species are found by 
simultaneously solving the differential equations, 
the conservation of atomic species relations, the 
fluid mechanics equations, and the electron energy 
equation. 

' Electron Temperature Equation 

As a consequence of the large ratio of atom or 
ion mass to electron mass, electrons transfer 
~ energy rapidly by collisions with other electrons 
but only slowly by elastic collisions with atoms or 
ions. The number of collisions necessary to pro- 
duce a Maxwellian velocity distribution between 
electrons and atoms is larger than the number of 
collisions necessary to produce a Maxwellian dis- 
tribution among electrons by a factor of the mass 
ratio. Therefore, throughout an ionizing shock 
wave, separate Maxwellian distributions describe 
the populations of electrons and heavy particles 



and thus, & different temperature la given to 
atoms and electrons in the same gas. 

The electron temperature is obtained from 
the solution to the electron energy equation. The 
change in electron energy is the sum of the energy 
exchanged by elastic and inelastic Collisions be* 
tween electrons and heavy particles, where in- 
elastic collisions are those which result In ioni- 
sation of the heavy particle 1 ! with resulting losses 
in election energy. A discussion of the electron 
energy equation for a nonequiHbrium gas is given 
by Appleton and Bray (Ref. 23). For a one- 
dimensional steady shock wave in a H^-He mix- 
ture, the resulting equation is 

~U e u) = T e ) 


Eq. (2) becomes 

1 + 3m“f a H°H (H) + ®He e He (He) 

( 3 ) 

where 

", ■ VV 4 "Jr'/”. * * R 4r>*V 

The elastic collision frequency can be written in 
terms of the elastic collision cross section using 
the approximate expression 

w ek " n k V e^ek 




where 


®He <*2 * 


R 2r + *4 


R 4r> 


(!) 


« molar concentration of species k, 


where t e » energy per unit volume of the electron, 
which is given by 

*. ■ }<•>**,• 

y = distance behind the shock wave, 

(e) = concentration of electrons, 

* mass of t pec lee k, 

v . * collision frequency for elastic colli- 
sions between electrons and species k, 

Qj' * ionization energy per mole of species k f 



and the (R^)'s the forward and reverse reac- 
tion rates for electron production. Equation (1) 
can be simplified using the steady- state approxi- 
mation, thereby setting the right-hand side of 
Eq. (1) equal to zero. This corresponds to the 
case where the energy is used by the electrons 
for ionization at the same rate that it is received 
from elastic collisions and has been discussed by 
Petschek and Byron (Ref. 1) and Hoffert and 
Lien (Ref. 5). The steady- state approximation 
results in the following relation for electron tem- 
perature: 


T * T 
e 


e H< R r R l/ R 3- R 3r> +e He< R 2- R 2r +R 4- R 4,) 


3m e V 


( 2 ) 


where the reaction rates are strong functions of 
T 0 due to the exponential dependence of the rate 
constants. An explicit equation for T» suitable 
for numerical solution can be obtained by approxi- 
mating the rate constants by linear functions of 
Ti for each integration step. Using the expres- 
sions 


and 


k H “ k l + k 3 “ *H T e + b H* 


k He " k 2 + k 4 “ a H« T e + b H«* 


v e * average electron velocity given by 

v e = (skT^mJ 1 / 2 , 

and 

Q . * the elastic collision cross section for 
e species k. 

The cross sections for electron-neutral collisions 
were obtained by fitting experimentally measured 
data with expressions of the form 

Q ek = A k exp *“ a E * " exp E *‘ 

The coefficients used are as follows: for He using 
data from Gould and Brown (Ref. 24), 

A = 7. 85 X 10‘ 16 , a = 0. 20, and b = 4.0; for H 2 
using data from Brode (Ref. 2 5), 

A = 2. 80 X 10*15, a - 0.3330, and b = 2.0; and 
for H using data from Brackmann, Fite, and 
Neynaber (Ref, 26), A = 3.95 X 10“15 f 
a = 0. 589, and b = a?. The Coulomb cross sec- 
tion given by the expression 

Q eH + = 1.95X 10* 6 T‘ 2 -°ln(l. 53 X 10 8 T 2 /n e )(cm 2 ) 
was used for H + and He + . 


C. Fluid Mechanics Equations 


The fluid mechanics of a one -dimensional 
shock wave with constant incident flow velocity is 
described, in a coordinate system fixed to the 
shock wave, by the continuity equation, 


Pl U l = P2 U 2* 

the momentum equation, 

Pj + PjU 2 = p 2 + p 2 u 2 , 
the energy equations, 

U 1 u 2 
h l + T * h 2 + 2 * 


(4) 


(5) 


( 6 ) 


A 


JPL Technical Memorandum 33-563 





I- 

% n 

‘fr- 





and the equations of state where the subscript (1) 
represents conditions in front of the shod' wave, 
subscript (2) represents conditions behind the 
shock wave, and p, u, and h are the mass density, 
fluid velocity, pressure, and enthalpy per unit 
mass respectively. For a multicomponent, chem- 
ically reacting gas with the electrons described by 
a separate temperature T e , the equations of state 
are 

Y. 

P = PRT= w 1 + n e RT e (7) 

and 

n_1 f 

h = 2 Y. (h7 + c T) + Y c T (8) 

flip. e p e 

1 * 1 r e 


where 'Wj is the molecular weight of species i and 
Y*^ = pi / p is the mass fraction of species i, h? is 
the hear of formation per unit mass of species i, 
and c is the specific heat per unit mass of spe- 
cies i . 1 


The value of the species concentrations are 
given by the relation 

£ (n./p) = R./p (9) 


where t is the particle time and Rj is the chemi- 
cal production rate of species i. It is not neces- 
sary to solve a rate equation for each species 
present In the system. An equation for the con- 
servation of atomic element may be written for 
each element present in the system. In terms of 
the frozen conditions, the conservation of atomic 
elements relation becomes 

2 a. .n./p = 2a..n.°/p\ (10) 

i ij i K i i 

where Ofy = number of atoms of element j per 
species x and (°) represents the frozen conditions 
in front of the shock wave. 


An explicit relation for the density ratio as a 
function of species mass fraction and T e is ob- 
tained by combining the conservation equations: 


_ -B - (b 2 - 4Ac) 1 '^ 2 

p 2 - 2A 


(ID 


where 


A 

B 

C 


r 


|m 2 ( 2 r - 1 ), 

-HI + Y M 2 - n e RT e / Pl ), 

r* + yM 2 /2 - r 1 (Y t - y«) h{ 
a l 

■ yYc T /a?, 

T e pe e' l 9 

n-1 

n_1 / 

R f Y/Wj 



M i* the Mach number, a U the incident efiaed of 
sound, v ifl the average specific heat, and the 
superscript (*) represents the froeen conditions 
in front of the shock wave, 

D, Continuum Intensity Calculation 

The continuum emission intensity during 
ionization relaxation is a function of electron con* 
centration, temperature, and the established 
radiative cross sections for hydroge.n Radiative 
processes that contribute to hydrogen continuum 
emission are radiative recombination 

H + + e - H + hv 

and radiative attachment of hydrogen 
H 4* e ->» + hv 

for temperatures below 1$,000*K, 

For an optically thin gas emission, intensity 
is given by the expression 

\ = B \(° r H n e ♦ a ir K H* n H n J l * 

where is the specific intensity in watt cm“ ' ster* 3 J, 
is the Planck function, or^ is the radiative re- 
combination cross section, <tjj- is the negative 
hydrogen ion absorption cross section, i = path 
length, is the equilibrium constant for H~, 

n e is the number density of electrons, and njj is 
the number density of hydrogen atoms. Empirical 
relations used for the cross sections at 
\ = 514. 5 rim we"e 

<r„ = (1.208 X 10~ 38 - 1.98 X 10‘ 43 T ) 

H e 

X exp (2. 86 5 X 10 4 /T e ) cm 2 

obtained from a fit to the absorption coefficients 
calculated by Stickford (Ref. 27) and 

<r„ = 2.21 X 10' 17 + 1.07 X 10‘ 21 T cm 2 
ri e 

obtained from a fit of the data presented by Griem 
(Ref. 28). The assumption has bien made that H“ 
is in local equilibrium with the electron concentra- 
tion during the course of the reaction. This is 
justified by the fact that the formation time for H~ 
is small compared with the ionization time. 

E. Results of the Ionization Rate Calculations 


The variation of species concentrations, tem- 
perature,, and emission intensity as a function of 
time have been obtained from numerical solutions 
to Eqs, (3), (9)» (10), and (11). The chemical 
rate equations were solved using a Runge-Kutta 
integration scheme with separate step size for the 
dissociation and ionization regions. In Fig. 3 
the variation in properties behind a shock wave 
are shown for a velocity of 1 5. 5 km/sec in a 0. 208 
- 0. 792 He mixture. The dissociation of H? it 
completed before appreciable ionization can take 
place. The increase in H concentration with time 
is due to the increase in total density during ion- 
ization, which offsets the removal of H atoms by 
ionization. The shape of the electron concentra- 
tion' rurve is similar to that obtained during argon 
ionization. Initially electrons are produced 
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primarily by Reactions 9 and 6. As the electron 
eimcentration and temperature increase, the rate 
of Reaction $ increases until it becomes the dom- 
inant electron producer. The arrow in Fig, 3 de- 
notes the point where the reaction rate fo*f Reac- 
tion 3 equals the rate for Reactions 5 and 6, 

-* rations 1, 2, a iv' 4 were found to be unimportant 
urn* r the conditions of this investigation but have 
been included for completeness. 

The electron temperature differs from the 
heavy particle temperature during the ionization 
relaxation, The initial drop in temperature is due 
to dissociation of and Is then followed by a 
slow increase in temperature during the early 
stages of ionisation. During the later stages of 
ionisation, the electron temperature asymptotic- 
ally approaches the heavy particle temperature. 

The deviation between the electron and heavy par- 
ade temperature was found to be proportional to 
the temperature behind the shock wave and there- 
fore the shock velocity. Results using the steady 
state treatment for the electron energy equation 
have been checked with exact numerical solutions 
of the electron energy equation, including a term 
for the energy of the electrons produced by atom- 


atom collisions. The steady state solution and the 
exact solution were found to converge within a 
time that was much shorter than the dissociation 
time. The atom- atom term and the initial value 
o l T e were found to have negligible effect on the 
solution a small distance behind the shock wave. 
This is consistent with estimates of the relative 
importance o' the atom-atom term using simpli- 
fied analysis. 


The significant difference found between T e 
and T during the ionization process is in dis- 
agreement with the results of Belozerov and 
Measures (Ref. 29). Their analysis indicated that 
the electron temperature in hydrogen ,closely fol- 
lowed the heavy particle temperature throughout 
the relaxation region and concluded that this was 
due to a large cross section for elastic H-e colli- 
sions. Their expression for the rate of energy 
transfer by elastic collisions contains an incor- 
rect integral of the cross section which overesti- 
mates the correct value by a factor of 40. Elec- 
tron temperature profiles in agreement with 
Belozerov and Measures can be produced using 
the expression 40 Q e n the present analysis. 
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IV. COMPARISON OF EXPERIMENTAL AND CALCULATED RESULTS 


Comparisons have been made between exper- 
imental measurements and computer calculations 
of the continuum intensity as a function of time. 
Important rate constants were varied in order to 
obtain a best fit with the experimental data and 
al*o to determine the sensitivity of the results to 
changes in each of the rate constants. From the 
agreement o' the calculation with the experimental 
data values for the atom- atom and atom -electron 
excitation rate, constants were obtained within the 
factor 01 uncertainty of the datti,. The reactions 
which govern the ionisation process are the atom- 
atom reactions (Reactions 5 and 6 } and the atom- 
electron excitation reaction (Reaction 3), while the 
direct ionization reactions were found to play a 
negligible role. In this investigation the assump- 
tion has been made that the collision cross sec- 
tions for electron excitation of hydrogen by H-H 
and F-He collisions are identical. Thus the rate 
constants *or Reactions 5 and 6 are identical, ex- 
cept for the reduced mass, and will be referred to 

as k AA* 

In Figs. 4 and 5 the effects of the atom-atom 
and electron-atom rate constants on the calculated 
continuum intensity are shown. The experimental 
emission intensity signal shown in Figs. 4 and 5 
has been normalized to the calculated equilibrium 
value. As can be expected, the atom -atom rate: 
constant affects only the early stage of the reaction. 
In Fig. 4 it is seen that an increase of the atom- 
atom rate constants has little effect on the maxi- 
mum slope of the intensity-time history but causes 
a steady decrease in the induction time (the inter- 
section of the maximum slope line with the abscis- 
sa). A best fit to the emission intensity signal is 
obtained for a value of the atom-atom rate constant , 
of n - 

^AA = 4 *° * \ - 

- ( 12 ) 

Using the values 0. 25 kAA and 4 kAA causes the 
calculated emission intensity to fall outside of the 
range of uncertainty of the experimental measure- 
ments. This is an indication of the sensitivity of 
the calculations to changes in the atom-atom rate 
constants. 

The effect on the calculations of variations of 
the electron^atom rate constant k^ is shown in 
Fig. The value for k^ was varied over the 
range of uncertainty in the electron beam meas- 
urement of electron-atom excitation cross sections 
while keeping the value of all other rate constants 
fixed. The maximum slope of the emission inten- 
sity and therefore the relaxation time was directly 


affected by tne value of ki, while the induction 
time was insensitive to changes in k^. The ex- 
pression for obtained from the cross-section 
measurements, 

k 3 = 7. 5X 10* 16 |^S.j 1 ^ 2 exp(-10.0AT e ) cm 3 sec' 1 , 

( 13 ) 

was found to give the best agreement with the 
emission intensity measurements. Due to the 
sensitivity of k^ to the induction time and k^ to 
the maximum slope, it was possible to unambigu- 
ously determine the expressions for each rate 
constant which best fit the experimental data. A 
comparison of a measured and calculated emis- 
sion intensity time history for another set of con- 
ditions, Uj = 13.7 km/sec and pj = 267 N/m^, 
is shown in Fig. 6 . 

The dependence of experimental and calculated 
relaxation data on shock velocity and initial pres- 
sure has been investigated. A characteristic time 
for ionization in the shock fixed coordinate sys- 
tem, T^ on , was defined as the product of the dens- 
ity ratio across the shock wave and tj a ^. In Fig. 

7, experimental values of (H2)2 T ion are shown 
a function of 1 /T 2 0 f where (H 2) 2 * s the H 2 CQI * 
centration immediately behind the shock wave ai i 
T 2 0 is the corresponding temperature. Also 
shown in Fig. 7 are calculated values of (H2)2 T ion 
for a range of values of kAA* The calculations 
are in good agreement with the experimental data 
using the expression for kAA g* ven in Eq. (12). 

The small dependence on p\ indicated by the cal- 
culated results was within the scatter of the data 
and could not be experimentally verified. Calcu- 
lated results using 0.25 kAA and 4. 0 kAA> shown __ 
in Fig. 7, bracket the experimental data and give 
upper and lower limits on the uncertainty in the 
determination of kAA* The error bars shown in 
Fig. 7 are indications of the experimental uncer- 
tainties in the determination of ionization time and^ 
the temperature. The largest uncertainty in T^ on 
occurs at the smaller values of I/T 2 n due to the 
approach of rise time limitations on the measure- 
ment of shock wave arrival time and radiation 
relaxation rate. For the slower shock velocities, 
the principal uncertainty in T£ cn is due to oscillo- 
gram reading errors cau ed by the decrease in 
signal- to -noise ratio. The uncertainties in the 
values for I/T 2 0 are caused by limitations on the 
accuracy of the shock velocity measurement and 
are largest at the lower shock velocities due to the 
increasing rate of change of I/T 2 0 with decreas- 
ing Uj. ’ ^ 
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V 


V. DISCUSSION AND CONCLUSIONS 


A model has been proposed for ionization 
processes governing shock structure in hydrogen- 
helium mixtures, and experimental data have been 
obtained which match the predictions of the model. 
The dissociation of hydrogen takes place in a 
region behind the shock wave that is small com- 
pared with ionization. The ionization process is 
governed by two separate regions: one dominated 
by the excitation of atomic hydrogen due to atom- 
atom collisions and the second due to hydrogen 
excitation by collisions with electrons. The elec- 
tron temperature was found to be significantly 
different from the heavy particle temperature dur- 
ing the early portion of the reaction but equaled 
the heavy particle temperature during the later 
part of the reaction. 

Data obtained on the relaxation time for the 
Hp signal is consistent with the postulated reaction 
model. The Hp relaxation times and the continuum 
relaxation times were found to be equal within the 
accuracy of the data. It was concluded from these 
results that the population of the n = 4 level of 
hydrogen follows the electron concentration. This 
is consistent with the reaction mechanism which 
assumes that all excited levels of hydrogen are in 
quasi- equilibrium with the electron concentration 
and the excitation of the first excited state is the 
rate limiting step. 

The agreement of the calculation with experi- 
mental emission measurements has produced esti- 
mates of the rate constants for the excitation of 
hydrogen to its first excited state by collision? 
with hydrogen and helium. The ratio of the cross 
section for hydrogen excitation by atom-atom 
collisions to the cross section for excitation by 
atom-electron collisions was found to be 
5. 34 X 10“2 b The cross section data are com- 
pared with data for other monatomic gases in 
Table 3, where values of S, the slope of the cross 
section curve at threshold, are shown for each 
species. It can be seen that while both the atom- 
atom and atom- electron excitation cross sections 
for hydrogen are considerably larger than those 
for the other gases, the ratio of atom-atom to 
atom-electron cross sections falls between the 
values for argon and krypton. The rate constant 
for atom-atom excitation of hydrogen was found to 
be on the order of 500 times larger^than^hat of 
argon, which is consistent with the finding the 
hydrogen-bearing impurities have a strong effect 
on argon ionization time measurements. 

The expression for k 3 that fits the data from 
this investigation is consistent with recent elec- 
tron beam cross section measurements of hydro- 
gen excitation. The rate constant expression for 
k^ given by Eq. (13) was obtained by combining 


the cross section data for excitation of the 2p 
level of McGowan, Williams, and Curley (Ref. 19) 
and the data for the excitation of the 2s level by 
Kauppila, Ott, and Fite (Ref. 20). This combined 
cross section is shown in Fig. 8 along with the 
combined cross section obtained using the 
McGowan et al. data and the 2s d^a of Lichten 
and Schultz (Ref. 31). Both curves for the com- 
bined cross section are seen to have the same 
initial slope. Error bars are used to show the 
experimental uncertainty in the measurements. 
Variations in the initial slope of the cross section 
analogous to the range of k 3 considered in Fig. 5 
are shown in Fig. 8. The upper cross section 
slope indicates the experimental uncertainty of the 
electron beam. The lower cross section slope 
was obtained using values from the earlier mea- 
surements of Fite, Stebbings, and Brackmann 
(Ref. 32), and Stebbings, Fite, Hummer, and 
Brackmann (Ref. 33). From Fig. 5, it can be 
seen that the sensitivity of the shock tube measure- 
ments for determining the excitation cross section 
is at least as good as the beam measurements. 

Thus shock tube data can be useful in obtaining 
cross section slopes at the threshold to supplement 
the data obtained from beam experiments. 

No results of ionization time measurements 
for hydrogen -helium mixtures have been pre- 
viously published. Studies of ionization rates in 
hydrogen or hydrogen-helium mixture using dia- 
phragmless electromagnetic shock-tubes have 
been made by Nakagawa and Wisler (Ref. 34), 
Wisler (Ref. 35), Fukuda, Okasaka, and Fujimoto 
(Ref. 36), and Belozerov and Measures (Ref. 29). 
Only Belozerov and Measures were able to pro- 
duce flow behind the shock wave that was thought 
to be free of driver gas mixing. Their results 
for the normalized ionization distance in hydrogen 
are shown in Fig. 9 along with results using the 
present analysis. The ionization distances calcu- 
lated using the expression for k^A from the pres- 
ent investigation are considerably larger than the 
measurements. However, calculated results using 
5. 0 k^A give good agreement with the results of 
Belozerov and Measures. This value of k/^A 
agrees with the present investigation if it is 
assumed that initial excitation of hydrogen occurs 
only by H-H collisions and that Reaction 6 is un- 
important. It was assumed in this study that the 
cross section for excitation was the same for both 
collision partners. At the present time there is 
no theoretical or experimental evidence to indicate 
that the excitation of hydrogen is much more prob- 
able by collisions with hydrogen atoms than by 
helium atoms. Subsequent studies of relaxation 
rates for a range of hydrogen- helium mixtures 
will be undertaken to determine the relative effi- 
ciency of these collision partners. 
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Table 1. Reaction scheme 


1. 

H + e 


H + + e + e 

z. 

He + e 


He + + e + e 

3. 

H + e 

z 

H* 4 e 

4. 

He + e 

■*- 

He* + e 

5. 

H + H 


H* + H 

6. 

H + He 

z 

H* + He 

7. 

H ? + He 

* 

H + H + He 

8. 

H 2 +H 2 

T 

H + H + H 2 

9. 

h 2 + h 

31 

H + H + H 

10* 

h 2 +h + 


H + H 4 H + 

11. 

H 2+ e 


H 4 H 4 e 
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Table 2. Rate constant expressions 


Reaction 


Rate constants m cm^ sec"^ 


Reference 


kj = 6.09 X 10 


-17 


/8kT X 1 / 2 

Uf) 


exp (-13. 6/kT e ) 


17 


k 2 = 3. 56 X 10 


-17 


C?f 


exp (-24. 58/kT e ) 


18 


= 7. 5 X lO- 1 *^) exp (-10/kT e ) 


19,20 


k. = 6.00 X 10 
4 


-17 




exp (-20/kT e ) 


21 


k 5 = 4.oxio- 17 (%3:y 


1/2 


exp (-10/kT) 


k 6 = 4.0 X 10- 17 (^1) 1/2 exp (-10/kT) 


k ? = 6. 93 X 10‘ 6 T _1 exp (103240/RT) 


11 


kg = 2. 5 k^ 


11 


k 9 = 20.0 k ? 


11 


10 


k 10 = k 9 


11 


k s lr 

K 11 k 9 


12 
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Table 3. Cross section slope at threshold for excitation by atom-atom and 
electron-atom collisions of monatomic gases 


Gas 

S AA’ Clr * 2 / ev 

S eA' 

2 / 

cm /ev 

S A A / S Ae 

Reference 

H 

4.0 X 10' 18 

7. 5 

-17 

X 10 1 

5. 34 X 

ID" 2 

Present in- 
vestigation 

Ar 

1. 2 X 10' 19 

4.9 

oo 

1 

o 

X 

2.45 X 

10- 2 

3 

Kr 

1.4 X 10' 19 

2. 32 

X 10‘ 18 

6. 0 X 

10" 2 

15 

Xe 

6.0 X 10' 21 

3. 80 

00 

1 

o 

X 

1. 58 X 

O 

1 

OJ 

15 

He 

2.0 X 10' 19 

2. 86 

X 10' 17 

7.0 X 

10- 3 

30 


4 

\ 


I 

I 



I 
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0 5 10 15 20 

TIME, /xsec 



TIME, /i$8c 

Fig, 3, Calculated variation of concentrations, 

and temperature with time behind a shock 
wave with shock velocity = 15. 5 km/sec 
and initial pressure = 133N/m^(1.0 
torr) 

14 


TRANSDUCER 

TIME ► 

Fig. 2. Emission intensity measur* ^nts for a 
shock velocity = 15.0 km/sec and initial 
pressure = 267 M/m 2 ( 2. 0 torr). The 
characteristic time for relaxation of the 
\ = 514. 5 nm signal is defined. 



0 12 3 4 

LABORATORY TIME, m mc 


Fig. 4. A measured continuum intensity-time 
history is compared with calculations 
using a range of values for the atom- 
atom rate constant. The solid line 
represents the experimental data. 
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LABORATORY TIME, n*c 

The effect of variations in 1^3 on the cal- 
culated continuum intensity-time history 


l ,n~12 




/// 

'//. 


//} 


!/T 2|0 *«0 5 , *K-' 


Fig. 7. Experimentally obtained normalised ion* 
Uatlon time compared with calculations^ 
A - Pl * 33.3 N/m 2 t o - p . - 133 N/m% 
O - pj * 267 N/m z . The dark curves 
represent calculations using the rate con* 
stant expressions shown in Table 2 for* 
from leftto right* pi ■ 33. 3* 133* and 
267 N/m 2 , The dashed curves indicate 
the effect of a factor of 4 increase and 
decrease to k AA , 


LABORATORY TIME, M *ec 

Fig. 6. Comparison of a calculated and measured 
intensity -time history for a run with 
slower shock velocity 


# 3 * 

- 0.3 


i/ 

\y 

r/ 


Fig. 8, Cross section for the excitation of hydro- 
gen by electron*atom collisions where 

- combined cross section 

from Refs, 19 and 20 


combined cross section 
from Refs. 19 and 30 

equivalent cross section 
for k 3 

equivalent cross section 
for 1.5 kj 

equivalent cross section 
for 0. 3 k 3 
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Uj, kftvfoc 


Fig. 9 


. Normalized ionization distance measure- 
manta of Belozerov and Measures, 

, compared with the re- 
sults of the present analysis using the 
rate constant expressions of Table 2 
_ _ • and 5. 0 k A A __ 
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